Dephasing in the adiabatic rapid passage in quantum dots: the role of phonon-assisted 
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We study the evolution of an exciton confined in a quantum dot adiabatically controlled by a 
frequency- swept (chirped) laser pulse in the presence of carrier-phonon coupling. We focus on the 
dynamics induced by a linearly polarized beam and analyze the decoherence due to phonon-assisted 
biexciton generation. We show that if the biexciton state is shifted down by a few meV, as is typically 
the case, the resulting decoherence is strong even at low temperatures. As a result, efficient state 
preparation is restricted to a small parameter area corresponding to low temperatures, positive 
chirps and moderate pulse areas. 
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I. INTRODUCTION 

The recently demonstrated^ high-fidelity preparation 
of a single exciton state in a self-assembled quantum 
dot (QD) by means of adiabatic evolution induced by 
a chirped laser pulse (referred to as adiabatic rapid pas- 
sage, ARP) opens new possibilities of charge control in 
QDs. In contrast to the more traditional Rabi flops^ as 
well as to methods based on voltage control^, ARP is 
much less sensitive to the details of the driving field. In 
particular, in the ideal case, the ARP technique ensures 
the full inversion of occupation between the empty dot 
and exciton states as soon as the pulse intensity reaches 
the threshold for the adiabatic passage. 

However, as QDs are embedded in a semiconduc- 
tor crystal matrix, carrier-phonon interactions are usu- 
ally found to considerably limit the fidelity of vari- 
ous optical control schemes. The phonon-induced de- 
phasing process is inevitable in such systems and leads 
to loss of information. These mechanisms have been 
investigated theoretically 6 8 and have been confirmed 
experimentally 9,10 . As we have recently shownii, in the 
case of the ARP the role of phonon coupling in the 
state preparation is essential. As the adiabatic evolu- 
tion follows one of the spectral branches correspond- 
ing to dressed exciton states, carrier phonon coupling 
leads to transitions to the other branch. Depending on 
whether the upper or lower branch is followed (which, in 
turn, depends on the direction of the chirp), the tran- 
sition is assisted by phonon emission or absorption, re- 
spectively. As a consequence, at low temperatures the 
phonon-induced decoherence leads to a strong asymme- 
try in the final occupation of the exciton state depending 
on the sign of the chirpii. While the analysis performed 
in Ref. [ll| was based on a two- level model, corresponding 
to the excitation with a circularly polarized laser pulse, 
in actual experiments linearly polarized beams have been 



used 1 . In this case, the selection rules allow the cou- 
pling to the biexciton state and phonon-assisted biexci- 
ton generation 12 may occur even at low temperatures, as 
the biexciton state is typically shifted to lower energies. 

In this paper, we study the evolution of the biexciton 
system driven by a linearly polarized chirped pulse in 
the regime of the ARP in the presence of the coupling to 
acoustic phonons. We show that including the biexciton 
state, which becomes relevant for linearly polarized laser 
pulses 1 , indeed opens a new decoherence path related to 
phonon-assisted transitions to the biexciton state. As 
a result, the fidelity of exciton state preparation drops 
down considerably in particular in most of the parame- 
ter areas where almost ideal evolution was predicted by 
the two-level model. Thus, high-fidelity control becomes 
restricted to the range of moderate pulse areas for posi- 
tive chirps and at low temperatures only, unless circular 
polarized pulses are used to suppress biexcitonic excita- 
tions. 

The paper is organized as follows. In Sec. HU we de- 
scribe the system and define the model used in our study. 
Sec. IIIII introduces our method of simulation of the open 
system dynamics. The results are presented in Sec. ITVl 
Finally, Sec. IVl concludes the paper. 



II. MODEL 

The biexciton system in the QD, restricted to the low- 
est bright states, is modeled as a four- level system with 
|0) representing the empty dot, being the two 

exciton states with different linear polarizations and \B) 
denoting the biexciton state. The corresponding Hamil- 
tonian is 

H x = E (\X)(X\ + \Y)(Y\) + (2E + E B )\B)(B\, 

where E$ is the exciton transition energy, Eb is the biex- 
citon shift, and we neglect the fine structure splitting 
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FIG. 1: (Color online) The adiabatic spectral branches of the 
optically driven QD system for G = 27T, to = 2 ps, \a\ = 
40 ps 2 and Eb — —2 meV: (a) positive chirp; (b) negative 
chirp. The evolution starting from the ground state follows 
the red line. The gray solid line represents the pulse envelope 
Q(t) which has been shifted upwards for clarity. 



which is irrelevant on the picosecond time scales consid- 
ered here. 

This charge system is driven by a chirped Gaussian 
pulse with the original envelope 



fto(*) 
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where is the original pulse area (before chirping) and 
To is the initial duration (pulse length). The central fre- 
quency uo is chosen at resonance with the single exciton 
transition. We assume that the linear chirp is effected by 
passing this pulse through a Gaussian chirp filter— with 
the chirp coefficient a. This results in the chirped pulse 
envelope 
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V / 27rr r 

with the chirped pulse length r 
the frequency chirp rate a = a/ (a 2 + t$ 
nian describing the coupling to the chirped pulse in the 
rotating wave approximation is 



~2r 2 
(a 2 /r 2 



+ T 2 ) 1 / 2 , and 
The Hamilto- 



Ml(t) 



(|0)(X| + \X)(B\) e ^ot+mi 2 /2 + h c 



Thus, for one linear polarization, the coupling is |0) O 
\X) <H> \B) and the other single-exciton state \Y) is 
decoupled (it is not coupled by phonon transitions, 
either—). Therefore, only a three- level model is needed 15 . 

The idea of the adiabatic passage^ is to drive the sys- 
tem state adiabatically along one of the spectral branches 
corresponding to the instantaneous eigenstates of the 
Hamiltonian Ho = Hx + H\ as . These instantaneous 
eigenstates are shown in Fig. [1] for the biexciton shift 
Eb = —2.0 meV. The branch followed by the system 
state is plotted with a thicker line. In Fig. DJa), we as- 
sume a positive chirp (a > 0), while in Fig. DJb) a neg- 
ative chirp is assumed. For pulse intensities above the 



ARP threshold, the central anticrossing at t = (be- 
tween the vacuum and exciton states) is passed adiabat- 
ically along a single spectral branch. On the contrary, 
the two other anticrossings away from t = (exciton- 
biexciton and vacuum biexciton) are met when the pulse 
is already very weak, hence they are very narrow (in the 
case shown here, the widths of both of these anticross- 
ings are below 10 /ieV, unresolved in Fig. [1]). At such 
narrow anticrossings, the sweep rate is very high com- 
pared to the distance between the two levels, which leads 
to a nearly purely non-adiabatic evolution: these anti- 
crossings are almost completely crossed by the evolving 
system state, that is, the system nearly completely jumps 
to the other branch (as shown by color coding of the lines 
in Fig.[T]). As a result, the charge configuration does not 
change here, as opposed to the central anticrossing. 

For the parameter range corresponding to the 
experiments 1,2 excitation of optical phonon modes is ex- 
cluded by the spectral properties of the system. In conse- 
quence, the coupling to acoustic phonons can be expected 
to dominate as a source of decoherence. The carrier- 
phonon interaction is described in terms of the usual in- 
dependent boson Hamiltonian, H int = S R, with 



and 



S=\X)(X\ + \Y)(Y\+2\B)(B\ 



R = ^2g k b k + h.c, 



where b k are phonon annihilation operators. We assume 
that the biexcitonic and excitonic wave functions can be 
factorized into products of single particle wave functions 
which are the same for both exciton and biexciton. This 
is well justified in the strong confinement limit where the 
wave functions are mostly determined by the confinement 
potentials. The coupling constants, accounting for the 
deformation potential interaction, can then be written 
as 



9k 



I hk 

2gvc 



[D e T e (k) - D h T h (k)}, 



where v is the normalization volume, g is the crystal den- 
sity, c is the speed of sound, D e ^) 18 the deformation 
potential constant for an electron (hole) and 



Fe(h)(k) = J d 3 r\^ e(h) (r)\ 



2 ^ikr 



are the form factors for electron (hole) wave function 
^e(h)( r )' Finally, the phonon reservoir is described by 
the Hamiltonian 



where uo k = ck are the phonon frequencies. 
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The system is conveniently described in the non- 
uniformly rotating frame defined by the unitary trans- 
formation e lA ^\ where 

A(t) = (u t + \at^ {\X){X\ + \Y)(Y\ + 2\B)(B\). 

The resulting Hamiltonian is H = H' x + H( as + H- m t+ -f/ph, 
with 

H' x = e iA ^Hxe~ iA{t) - KA(t) 

= hA(t) (\X)(X\ + \Y)(Y\) + (2hA(t) + Eb) \B)(B\, 

where A(t) = —at, and 

H{ as = e iA ^H las e~ iA ^ 

= ^(|0XX| + |XXB|) + h.c. 



III. METHOD OF SIMULATION 

To simplify the numerical simulation, we perform the 
unitary transformation defined by the operator 



(\X)(X\ + \Y)(Y\)W + \B)(B\W 2 , 



where 



W = exp 



h 



Upon this transformation and expanding to the leading 
order in the coupling constants, the Hamiltonian may be 
written as 



H = W#W T = fTi + H ph + V, 
where Hi = H' x + wH[^ s and 

V= ^l mx \ + \X)(B\)W^+h.c. 

Here 
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accounts for the phonon-induced renormalization of the 
pulse amplitude in the slow driving limit 7,10 . We assume 
that the central frequency of the driving pulse is corrected 
for the phonon-induced energy shifts. The evolution of 
the reduced density matrix of the charge subsystem is 
then found by numerically solving the lowest-order time- 
convolutionless (TCL) evolution equatio n 17 ! 18 



-f 

Jo 



drTr ph [V{t),[V{r),p{t)®p^]\, 



where V(t) is the interaction Hamiltonian V in the in- 
teraction picture with respect to H\ + iJ p h, p P h is the 
phonon density matrix at the thermal equilibrium, and 
Tr p h denotes the partial trace with respect to phonon 
states. 
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FIG. 2: (Color online) The final occupation of the single- 
exciton state (color coded) as a function of the original pulse 
area G and the chirp a in the absence of coupling to phonons: 
(a) for E B = -1.0 meV (b) for E B - -2.0 meV. 
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FIG. 3: (Color online) The final occupation of the single- 
exciton state (color coded) as a function of the original pulse 
area G and the chirp a with phonon effects included: (a,c,e) 
for E B = -2.0 meV (b,d,f) for E B = -8.0 meV. 



IV. RESULTS 

With the simulation method outlined in Sec. IIII| we 
study the exciton dynamics focusing on the possibility of 
phonon-assisted transitions to the biexciton state, when 
the system is excited with a linearly polarized laser pulse. 
In our simulations, we assume g = 5360 kg/m 3 , c = 
5110 m/s, D e = 7 eV and D h = -3.5 eV. The electron 
wave functions are assumed to be Gaussians with the 
extension l e = 4.5 nm in the xy plane and l GyZ = 1 nm 
along the growth direction. In the case of holes, we take 
lh = 0.87/ e and lh, z — 0-87/ e ,z respectively. The pulse 
length (before chirping) is taken to be To = 2 ps. 

First, we need to assess the degree of perturbation to 
the adiabatic evolution which stems from the additional 
level crossings with the biexciton (see Fig. [1]) and is not 
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FIG. 4: (Color online) The final occupation of the single- 
exciton state as a function of the original pulse area G at 
two different temperatures: (a,b) for Eb = —2.0 meV; (c,d) 
for Eb = —8.0 meV. Blue dotted line shows the results with 
chirp a = — 20 ps 2 , the red solid line represents a = 0, while 
the green dashed line corresponds to a — 20 ps 2 . The final 
occupation of the biexciton state for Eb — — 2.0 meV is shown 
in (e,f). 
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FIG. 5: (Color online) The total number of excitons (color 
coded) as a function of the original pulse area G and the chirp 
a for Eb = -2.0 meV (a) at T = 1 K and (b) at T = 20 K. 



related to phonons. To this end, in Fig. [2] we show the 
single-exciton occupation without carrier-phonon cou- 
pling for Eb = —1.0 meV and Eb = —2.0 meV. The 
results show a reduction of the final occupation of the 
single exciton state at higher pulse areas for a very small 
biexciton shift (Eb = — 1.0 meV) and positive chirps. We 
have found that under these conditions, a biexciton oc- 
cupation on the order of 0.1 occurs (not shown). On the 
other hand, for the still small value of Eb = — 2.0 meV, 
the evolution is nearly unperturbed by the presence of the 



biexciton state. Thus we conclude that for \Eb \ ^2 meV 
no direct population of the biexciton occurs with the 
pulses studied here, as the driving of the exciton to biex- 
citon transition is too off-resonant. 

In the next step, we include phonons and consider two 
cases: a) the biexciton state is shifted by a typical energy 
(Eb — — 2.0 meV) below the single-exciton state, b) for 
comparison, the biexciton is shifted by a very large (and 
rarely realized) energy (Eb = —8.0 meV). In the first 
case, one expects that phonon-induced transitions to the 
biexciton state are efficient. In the second case, due to 
the large energy separation between the states (beyond 
the cut-off frequency of the carrier-phonon coupling), the 
biexciton generation is expected to be negligible and the 
system can effectively be treated as a two-level system. 
The simulated evolution of the single exciton occupa- 
tion for both cases is shown in Fig. [3j The left panels 
(Fig. [U[a,c,e)) contain results for Eb = —2.0 meV, and 
the right panels (Fig. EP>,d,f)) for E B = -8.0 meV. The 
final occupation is calculated as a function of the pulse 
area 6 and the chirp rate a at three temperatures (T = 1, 
20, and 40 K). 

Let us start the discussion for the case of the large biex- 
citon shift Eb = —8.0 meV and low temperature T = 1 K 
shown in Fig.[3jb). As expected, the results from the two- 
level model 11 are retrieved and a strong asymmetry in the 
final occupation with respect to the chirp is seen. In the 
two-level model, for negative chirp the evolution follows 
the upper branch and phonon emission leads to damping, 
while for positive chirps the lower branch is followed and 
since phonon absorption is almost absent at 1 K a full ex- 
citon occupation is reached. In contrast, in Fig.[3fa), for 
the small biexciton shift Eb = — 2.0 meV, the influence of 
the biexciton state can be clearly seen. In particular, for 
positive chirps a decoherence at high pulse areas is vis- 
ible due to the opening of the new decoherence channel 
consisting in transitions to the biexciton state (phonon- 
assisted biexciton generation). As can be seen in Fig. [TJ 
for a negative biexciton shift (positive biexciton binding 
energy), this process corresponds to phonon emission no 
matter which adiabatic branch is followed, hence it is ef- 
fective irrespective of the sign of the chirp even at low 
temperatures. Also for negative chirps the damping has 
become more effective. Thus, only for moderate pulse 
areas around = 2ir a high-fidelity preparation of the 
single exciton state is possible. 

When considering higher temperatures, as seen in 
Fig. Efc-f), phonon absorption becomes more likely and 
now also for positive chirps decoherence due to transi- 
tions between the main branches is stronger. Hence, con- 
siderable phonon-induced perturbation is seen even in the 
case of strongly detuned biexciton (Fig. |2fd,f)). In con- 
sequence, the additional effect of the biexciton shift is 
relatively weak. 

For a more quantitative analysis we show the final oc- 
cupation as a function of pulse area for different chirps 
a = -20, 0, and 20 ps 2 in Fig.^a-d) for E B = -2.0 meV 
(Fig.m^b)) and for E B = -8.0 meV (Fig. E(c,d)). The 
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most remarkable difference for the two biexciton shifts 
is again seen for positive chirp a = 20 ps 2 at T = 1 K. 
While in Fig. [4f a) , for the small biexciton shift, a strong 
damping of the exciton occupation is found, in Fig.HJc), 
for the large biexciton shift, an exciton occupation of 
one is reached for all pulse areas above the adiabatic 
threshold. Also the decoherence for the negative chirp 
a = — 20 ps 2 and damping of the Rabi oscillations at 
zero chirp are stronger due to the phonon-assisted biex- 
citon generation. For higher temperature at T = 40 K in 
Fig. |4^b,d) the difference between small and large biex- 
citon shift is less pronounced. 

The strong decoherence for the small biexciton shift 
Eb = —2.0 meV is accompanied by an occupation of 
the biexciton state as shown in Fig.|4je,f) confirming our 
previous statements. 

Since in some experiments 2 a photocurrent detection 
scheme is used, in which the detection signal is propor- 
tional to the total number of photocreated excitons, we 
have calculated also the total number of excitons in the 
system, iV tota i = A/x+2A/"b, where Nx and TVb are the oc- 
cupation of the single- and biexciton state respectively. 
The results are shown in Fig. [5] for Eb — —2 meV at 
T = 1 K (Fig. GJa)) and T = 20 K (Fig. GJb)). One can 
see that at a certain level of phonon-assisted biexciton 
generation the result exceeds 1, which has no sense if in- 
terpreted in terms of the two-level model. Moreover, it 
is clear that even below this limit, the biexciton contri- 
bution leads to an overestimated occupation of the single 
exciton state. 



V. CONCLUSIONS 

In summary, we have studied the open system evolu- 
tion of a QD adiabatically controlled by a chirped linearly 
polarized laser pulse. We have shown that the presence 
of a biexciton state, which becomes optically coupled for 
linearly polarized laser pulses \ opens new decoherence 
paths related to phonon-assisted biexciton generation. If 
the biexciton state is shifted down by 1 or 2 meV it be- 
comes accessible via phonon emission processes and the 
resulting decoherence is strong even at low temperatures. 
This reduces the achieved occupation of the single exci- 
ton state in the parameters range where the chirped pulse 
control was very efficient in a model without the biexci- 
ton state. Since the phonon-assisted biexciton generation 
increases with the pulse area moderate pulse areas are 
optimal for the state preparation. On the other hand, 
if the biexciton shift is larger, the two- level evolution is 
essentially recovered. 
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